Relativistic unitary coupled cluster theory and applications 
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We present the first formulation and application of relativistic unitary coupled cluster theory to 
atomic properties. The remarkable features of this theory are highlighted, and it is used to calculate 
the lifetimes of 5 2 D3/2 and 6 2 P3/2 states of Ba + and Pb + respectively. The results clearly suggest 
that it is very well suited for accurate ab initio calculations of properties of heavy atomic systems. 
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There have been a number of attempts to modify 
coupled-cluster (CC) theory 0, despite its spectacular 
success in elucidating theproperties of a wide range of 
many-body systems [H 0, 0, 0, Q • One interesting case in 
point is unitary coupled-cluster (UCC) theory which was 
first proposed by Kutzelnigg (| . In this theory, the effec- 
tive Hamiltonian is Hcrmitian by construction and the 
energy which is the expectation value of this operator 
in the reference state is unlike in CC theory, an upper 
bound to the ground state energy Another attrac- 
tive feature of this theory which we shall discuss later 
is that at a given level of approximation it incorporates 
certain higher order excitations that are not present in 
CC theory. Furthermore, it is well suited for the cal- 
culation of properties where core relaxation effects are 
important |3> Ig • I n spite of the aforementioned advan- 
tages, there have been relatively few studies based on this 
method [jj ^| . This work is the first relativistic formu- 
lation of unitary coupled cluster (UCC) theory and also 
the first application of this theory to atomic properties. 
In this letter, we first present the formal aspects of rel- 
ativistic UCC theory and then apply it to calculate the 
lifetimes of the 5 2 D 3 / 2 and & 2 P^/2 states of Ba + and Pb + 
respectively; which depend strongly on both relativistic 
and correlation effects. The comparison of the results of 
these calculations with accurate experimental data would 
constitute an important test of this theory. 

The exact wave function for a closed shell state in CC 
theory is obtained by the action of the operator exp(T) 
on the reference state |$). However, in UCC theory jlCj . 
it is written as 



*) = exp(cr) |$) 



(1) 



where a = T - Tt; and T and are the excitation 
and deexcitation operators respectively, a is clearly anti- 
Hermitian, since <j* = — a. Using this unitary ansatz for 
the correlated wave function, the relativistic UCC equa- 



tion in the Dirac-Coulomb approximation can be written 

as 

exp(o- t )iZ'exp(o-)|$) =J5|$), (2) 
where H is the Dirac-Coulomb Hamiltonian 



H = Y,ca t - Pl + ((3-l)mc 2 + V N + Y,— ■ (3) 



Using the normal ordered Hamiltonian, Eq.|J2J can be 
rewritten as 



exp(cr t ) J ff A r exp(cr) |$) = AE |$) , 



(4) 



where the normal ordered Hamiltonian is defined as 
H N = H - ($| H |$) and AE = E - ($| H |$>. . 

The choice of the operator a makes the effective Hamil- 
tonian Hn = cxp(— u)Hn exp(cr) Hermitian. 

The effective Hamiltonian is expressed by the Haus- 
dorff expansion in CC theory as 

H N = exp(-T)H N exp(T) 

= H N + [H N ,T] + ^[[H N ,T],T] 

+ ^[[[H N ,T],T),T] 

+ ^[[[[H N ,T],T],T},T}. (5) 

In UCC, the operator T is replaced by a = T — in the 
above equation and this results in Hm being expressed 
in terms of a non-terminating series of commutators. For 
practical reasons, one truncates the series after some fi- 
nite order. Truncation at the n-th order commutator 
leads to the nomenclature UCC(?t.). 

Using UCC( 3) approximation and without modifying 
the last term of Eq. we write Eq. J3} as 
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H N + H N T + - ( .//,, TT : -27T//. T 



=)4 



H N TTT + 3T^H N T 



3T^H n TT\ + —HffTTTT 



|$) = AE\$) 
(6) 



By projecting single, double and higher order excited 
determinant states on Eq.©, we get the cluster ampli- 
tude equations. The approximation which includes only 
single and double excitations/deexcitations in the UCC 
wavefunction is known as UCCSD. After careful analysis, 

one finds that there are new terms arising from T^HnT, 



TWH n T and T^H N TT in the UCCSD equations. The 
first term will give rise to some extra diagrams which cor- 
respond to double and triple excitations and the last two 
terms to triple and quadrupole excitations. Some typical 
diagrams which represent triple and quadrupole excita- 
tions are given in figures and \T]p respectively. The 
UCCSD approximation in addition to single and dou- 
ble excitations, also includes some triple and quadrupole 
excitations to all orders in the residual Coulomb inter- 
action in a more elegant and simpler manner compared 
to CC theory. It would have been computationally pro- 
hibitive to include triple and quadrupole excitations for 
heavy atoms in the framework of normal CC theory. In- 
deed one of the principal advantages of UCC theory is its 
ability to subsume higher levels of excitations than CC 
theory at the same level of approximation. 

To calculate the ground state of the system we 
first compute the cluster amplitude of the closed shell 
states(Ba ++ and Pb ++ ) by using Eq.© and then use 
the open shell coupled cluster method (OSCC) 11]. The 
exact wave function can then be written as 



N+l 



)=exp( ( r){l + ^}K +1 ) 



where | < l , ^ r+1 ) is the Dirac-Fock reference state which we 
get after adding an electron to the fcth virtual orbital and 
Sk is the corresponding excitation operator. We obtain 
a set of equations 



N+X \Hn(1 



H, 



cir 



(8) 



and 



k n+i 



H N (l + S k )\^ +1 ) 



N+l 



s k \* 



N+l 



)H eS . 



(9) 



The Eq. © is non- linear in S k because H c s is itself a func- 



tion of S k where S k = J2ak s a a W + Ea&fcr s ab a k a t a ba a ; 

and s^l are the single and double excitation cluster 
amplitudes for the valence electrons. The labels a, b and 
k,r refer to the core and virtual orbitals respectively. 
Hence, these equations have to be solved self-consistently 
to determine the S k amplitudes. We have included the 
triple excitations in our calculations in an approximate 



way. The amplitudes corresponding to these excitations 
are of the form 



opqr _ 

°abk — 



VT 2 + VS 2 



Sb + £k — £p — £q 



(10) 



where S^? k are the amplitudes corresponding to the si- 
multaneous excitation of orbitals a, b, k to p, q, r respec- 
tively and VT and VS are the correlated composites in- 
volving V and T, and V and S respectively, e's are the 
orbital energies. The above amplitudes are added ap- 
propriately in the singles and doubles open shell clus- 
ter amplitude equations and they are then solved self- 
consistently. We therefore obtain solutions of 5*i and S 2 
amplitudes taking into consideration the partial effect of 
the triple excitations. We had referred to this approxi- 
mation earlier as CCSD(T). 

As we have seen in Eq.(jHJ, Hat is defined as Hn = 
exp(— (t)Hn exp(cr). From Eq.© it is clear that the the 
expansion of H jv corresponding to UCC theory can be 
expressed as 



H 



N 



(Hn) 



JVjCCSD 



(Hn) e 



(ii) 



Following much the same argument as in the closed shell 
case, it can be shown that the contributions from the ex- 



tra part arising in Ea. (|ll|) contains terms like T^HjyT, 

TtTtiJ^T and T^TTnTT. These terms will affect the 
open-shell amplitude determining equations. In the UCC 
theory for the closed shell, we have seen that although 
there are some extra terms arising from the expansion of 
Hn, the contributing diagrams corresponding to those 
terms are not always new. Hence we have taken into ac- 
count those diagrams which are not present in the CCSD 
approximation. It can be shown that the extra terms 
arising from the inclusion on the operator a = T — T** 
give rise to some new terms but the corresponding dia- 
grams have already been taken into account through the 
inclusion of partial triple amplitudes. In analogy with 
CCSD(T), we refer to our approach as UCCSD (T). 

The normalized transition matrix element (i — ► /) 
due to an operator O is given by 
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We consider next the calculations of the lifetimes of the 
5 2 D 3 /2 and Q 2 Pz/2 states of Ba + and Pb + respectively. 
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FIG. 1: Typical triples quadruples diagrams arising from UCC(3) 
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FIG. 2: Decay scheme of the low lying states of Ba + and Pb + 



Both of these states decay to the ground states via Ml 
and E2 transitions as shown in figure [21 The transition 

probabilities A (Aft — A t >f in s _1 ) for Ml and E2 

transitions are expressed as 1 1 121 ] 



AMI 



2.6973 x 10 1 



qMI 
~ a fi 



and 



E2 _ 1.1199 x 10 18 oE2 



(12) 



(13) 



respectively, where Ji is the degeneracy of the initial state 
and A (in A) is the wavelength corresponding to the tran- 
sition i — > /. In Eas. (|12fl and (|13() the line-strength 
Sfi is |Ml/i| and |J52/i| respectively and Ml/j and 
E2 fi are the corresponding one electron reduced matrix 
elements of magnetic dipole and electric quadrupole tran- 
sitions [lflj . 

The net probability for a given transition which allows 
two different channels is given by |l^ 



TABLE I: Description of the basis functions used in the UCC 
calculation of Ba + . 





Sl/2 Pl/2 PS/2 


d-i/2 


^5/2 


h/2 


77/2 57/2 ff9/2 


Analytical 


5 6 6 


5 


5 


8 


8 6 6 


Numerical 


8 7 7 


5 


5 


2 


2 


Total 


13 13 13 


10 


10 


10 


10 6 6 



TABLE II: Excitation energies (in cm 1 ) and lifetime (in 
s) of 5 2 D s /2 state of Ba + corresponding to the transition 



5 Z D. 



3/2 



/2 





Excitation energy 


r 


UCCSD(T) 


4789 


81.01 


CCSD(T) yj 


4809 


87.06 


Dzuba et al [15] 


4411 


81.5 


Guet et al [16] 


4688 


83.7 


Experiment ^JJ 


4874 


79±4.6 



A = A 



Ml 



A 



E2 



(14) 



and the corresponding lifetime, which is the inverse of 
the transition probability is expressed as 



1 

TMl 



1 

TE2 



(15) 



The results of the calculations of the lifetimes of the 
5 2 D 3 / 2 state of Ba + and 6 2 P 3 / 2 state of Pb + are given 
in tables [H] and IIIII respectively. These calculations are 



very challenging as they involve accurate determinations 
of the line strengths of the Ml and E2 transitions and 
more critically the third and fifth powers of the excitation 
energies. 

The lifetime of the 5 2 D 3 / 2 state of Ba + is calculated 
using hybrid orbitals consisting of single particle orbitals 
which are partly numerical and partly analytical. Such 
an approach is described in detail by Majumder et al 
[Tsj . The analytical orbitals are Gaussian type orbitals 
(GTOs). The details of this basis are given in table HI 
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TABLE III: Excitation energies (in cm 1 ) and lifetime (in s) 
of the 6 2 P 3/2 state of Pb + 





Excitation energy 


T 


Dirac-Fock 


13612 


0.0415 


CCSD(T)[13] 


13710 


0.0425 


UCCSD(T) 


13719 


0.0413 


Exp [20] 


14085 


.0412(7) 



Excitations from all the core orbitals are included for 
this ion as well as Pb + . Since the calculation of transi- 
tion probabilities/lifetimes are extremely sensitive to the 
excitation energies. Gopakumar et al |14| have shown 
that the hybrid basis can serve this purpose very well. 
Our UCCSD(T) calculations yield ionization potentials 
(IPs) to an accuracy of about 0.1%, and 0.07% for the 
6 2 <Si/2 and 5 2 D 3 / 2 states of Ba + respectively. The cor- 
responding IP (in cm -1 ) values are 80544 (80687) and 
75755 (75813). The numbers in the parentheses are the 
corresponding experimental values . 

The values of the calculated and experimental lifetimes 
are given in table |H] Dzuba et al\H^, Guet et al 
and Gopakumar et al have not considered the Ml 
channel in their calculations. Gopakumar et al have 
used CCSD(T), but Guet et al and Dzuba et al have 
used different variants of many-body perturbation the- 
ory (MBPT) with certain semi-empirical features to cal- 
culate the transition amplitudes. In addition, the latter 
two have used experimental values of the excitation en- 
ergies to calculate the transition probabilities/lifetimes. 
Our calculation based on UCCSD(T) is purely ab ini- 
tio and takes into account the contribution of the Ml 
transition. The transition amplitude in this approach 
includes some terms in addition to those that appear in 



CCSD(T); and they arise due to the presence of the deex- 
citation operator in the UCC wavefunction. The result of 
the lifetime calculation considering only the E2 channel 
is 85.567. Inclusion of the Ml transition improves the to- 
tal value (81.01 s) which is within the experimental error 
bar (79±4.6 s). 

The UCCSD(T) calculations of lifetime for the 6 2 P 3/2 
state of Pb + are carried out using the same basis func- 
tions used by Sahoo et al 0|- The leading correlation 
contribution to the Ml and E2 channel come from the 
core polarization and pair correlation effects. It is evident 
from Table Hill that the accuracy of the calculation of ex- 
citation energies using UCCSD(T) is even better than 
that obtained from CCSD(T). This result and the im- 
provements in the Ml and E2 transition amplitudes due 
to the UCC formulation coupled together give a value 
of the lifetime of the 6 2 P 3 / 2 state of Pb + that is clearly 
in better agreement with the corresponding CCSD(T) 
calculation. This value is within the limits of the exper- 
imental error [2£]], but this is not the case for CCSD(T). 

In summary, we have applied the relativistic UCC the- 
ory for the first time to atomic properties. The results of 
our UCCSD(T) calculations of the lifetimes of the 5 2 Z) 3 /2 
and 6 2 P 3 / 2 states of Ba + and Pb + respectively are in 
very good agreement with experiment and superior to 
those of all previous calculations. It indeed appears that 
this theory is capable of yielding high precision results for 
a wide range of properties of heavy atomic systems in- 
cluding violation of parity and time-reversal symmetries. 
In addition, it would be worthwhile to explore its feasi- 
bility for the accurate determination of the properties of 
other many-body systems, particularly when it becomes 
necessary to go beyond the usual coupled cluster theory. 
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